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A density wave in quark matter is discussed at finite temperature, which occurs along with the 
chiral condensation, and is described by a dual standing wave in scalar and pseudo-scalar condensates 
on the chiral circle. The mechanism is quite similar to that for the spin density wave suggested by 
Overhauser, and entirely reflects many-body effects. It is found within a mean-field approximation 
for NJL model that the chiral condensed phase with the density wave develops at a high-density 
region just outside the usual chiral-transition line in phase diagram. A magnetic property of the 
density wave is also elucidated. 

I. INTRODUCTION 

In a recent decade the condensed matter physics of QCD has been an exciting area in nuclear physics. In 
particular, its phase structure at finite-density and relatively low-temperature region is studied actively, since 
it was suggested that the color superconductivity (CSC) involved observable consequence of quark matter due 
to the large magnitude of its gap energy over a few hundred MeV P10,OJ- Although QCD at finite densities 
has important implications for physics in other fields, e.g., the core of the compact stars and their evolution Q, 
it remains in not understood adequately. 

Whereas the quark Cooper-pair (p-p) condensations attracts much interest, particle-hole iv-h) condensations, 
which are related to the chiral condensation 0, i or ferromagnetism in quark matter B & 0, have also 
been studied, and their interplay with CSC has been discussed. Since the Cooper instability on the Fermi 
surface occurs for arbitrarily weak interaction, the p-p condensation should dominate at asymptotically free 
high-density limit. While on the other hand, there exists a critical strength of the interaction in the p-h channel 
for its condensation. 

At moderate densities, however, where the interaction is strong enough, p-p and p-h condensations are com- 
petitive, and various types of the p-h condensations are proposed j3,|J in which the p-h pairs in scalar or tensor 
channels have a finite total momentum indicating standing waves. Instability for the density wave in quark 
matter was first discussed by Deryagin et al. [lOj at asymptotically high densities where the interaction was 
very weak, and they concluded that the density-wave instability prevailed over the Cooper's one in the large N c 
(the number of colors) limit due to the dynamical suppression of colored p-p pairings. 

In general, the density waves are favored in one-dimensional (1-D) systems, and have the wavenumber Q = 2kp 
according to the Peierls instability [ill H^ . e.g., charge density waves in quasi- 1-D metals ^3|. The essence of 
its mechanism is the nesting of Fermi surfaces and the level-repulsion (crossing) of single particle spectra due 
to the p-h interaction with the finite total wavenumber. Thus the low dimensionality has a essential role to 
make the density-wave states stable. In higher dimensional systems, however, the transition occurs provided 
the interaction of a corresponding (p-h) channel is strong enough. For the 3-D electron gas, it was shown by 
Overhauser [LH 1 1 ~J| that the paramagnetic state was unstable with respect to formation of a static spin density 
wave (SDW), in which spectra of up- and down-spin states deform to bring about a level-crossing due to the 
Fock exchange interactions, while the wavenumber does not precisely coincide with 2kp because of incomplete 
nesting in the higher dimension. 

In the recent paper , we suggested a density wave in quark matter at moderate densities in analogy with 
SDW mentioned above. It occurs along with the chiral condensation and is represented by a dual standing wave 
in scalar and pseudo-scalar condensates (we have called it 'dual chiral-density wave'-DCDW). The DCDW has 
different features in comparison with the previously discussed chiral-density waves 0, 0, IT(i| , and emerges at a 
moderate density region pb/po — 3 — 6 (where po = 0.16fm -3 , the normal nuclear matter density). 

In this paper we would like to further discuss DCDW and figure out the mechanism in detail. We also present 
a phase diagram on the density-temperature plane. In Sec. II we start by introducing the order parameters for 
the dual chiral-density wave and show the nature of the ground state on the analogy of the spin density wave: 
a level-crossing of single-particle spectrum in the left- and right-handed quarks. Sec. Ill is devoted to concrete 
calculations by use of Nambu-Jona-Lasinio (NJL) model [l7| to present a phase diagram at finite density and 
temperature, in the case of 2 flavors and 3 colors. In Sec. IV we summarize and give some comments on 
outlooks. 
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FIG. 1: Two branches of the single-particle spectrum are plotted along p z axis at p x — p y — where q/ jz. Dashed 
(Solid) lines correspond to M = (M ^ 0). 



II. NATURE OF DUAL CHIRAL-DENSITY WAVE 

In the pioneering studies of chiral density waves H U ^| > a spatial modulation in the chiral condensation 
was considered; the scalar condensation with a wavenumber vector q occurs, {ipip) <x cos(q • r). In this section, 
we consider a directional modulation with respect to the chiral rotation, since the directional excitation modes 
should be lower than the radial ones in the spontaneously symmetry-broken phases. We propose the dual 
chiral-density wave (DCDW) in scalar and pseudo-scalar condensation, 

(tyij)) = Acos(q-r), 
(V^7 5 i/>) = Asin(q-r), (1) 

where the amplitude A corresponds to the magnitude of the chiral condensation; {i)ijj) 2 + (ipi'ys'ip) 2 = A 2 . 

We give the single-particle spectrum, in the presence of the density wave as an external field. The 
Lagrangian density of the system in the chiral limit then reads 

C = VJ(r)[i^ + 2GA{cos(q-r) + i75sin(q-r)}]^(r), (2) 

where 2G is a coupling constant between the quark and DCDW. In the chiral representation, it becomes clear 
how the density wave affects the quark fields: the Lagrangian shows that the density wave connects left- and 
right-handed particles in pairs with the total momentum q, 

A\r - [ d3p ( ^(P-g/2)V ( ^(p-g/2)* —M \(i> L (p-q/2)\ () 
J (2tt)3 {iP R (p + q/2) J { —M a^p + q/2Y ) \ ^ R ( P + q/2) ) ' [6) 

where q^ — (0, q), and M(= —2GA) corresponds to the dynamical mass if the density wave is generated by 
quark interactions as the mean-field. The single-particle (quasi-particle) spectrum is obtained from the poles of 
the propagator: 

det [a^p - g/2)"oy,(p + g/2)" - M 2 ] = 0. (4) 
From the above equation we found the spectrum: positive and negative energies, E±(p) and — E±(jp), 



E±(p) = V^o(p) 2 + |q| 2 /4 ± V(P ■ q) 2 + M 2 |q| 2 , E a (p) = (M 2 + Ipl 2 ) 1 / 2 . (5) 

Because of the finite q and A, the spectrum is deformed and split into two branches denoted by E±(p), as 
shown in Fig. ^ where the direction of q is taken parallel to the z axis. The two branches exhibit a level-crossing 
between the left- and right-handed particles. 

For q = (0, 0, q), the corresponding eigen spinors are given, in the chiral representation, by 



^(p-q/2) 
^(p + q/2) 

V±(p) = it±(p)|E ± ( P )^-B±(p) for negative energy- E±(p), (7) 



u±(p) = N± ( fi), , ,J ) for positive energy E±(p), (6) 
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where Af± is a normalization factor for Uj s1 u± = v±v± = 1, and 

^ ±} (p + q/2) 
^(p-q/2) 

Here it should be noted that the original spinor in eq. @, ip(r), is represented as a plane wave expansion for 
the eigen spinors in eqs. ©-10: 

^(r) = e -^q-r/2 J- J J^L {a,(p)u a (p) + b s (p)v a (p)} e~ ip ' r = e+i***l* J J^L^e*" , (10) 

where a±(p) {b±(p)) is annihilation operator for the positive (negative) quasi-particle state. The factor, 
exp(— ijsq ■ r/2), comes from the momentum shifts in the left- and right-handed particles in eq. © 

due to the presence of DCDW, and thus reflects a nature of the spatially chiral-rotated ground state. The fac- 
tor corresponds also to a kind of Weinberg transformation, which changes the system from spatially modulated 
one to a uniform one 0, 0] . 

In the previous section, we mentioned that in general, density waves should be favored in 1-D systems. If we 
assume a quasi- 1-D system along the direction of the z axis, suppressing the radial (x-y) degrees of freedom, a 
gap (~ 2M) opens just above the Fermi surface, provided that q is taken to be 2k F (kF' the Fermi momentum 
for free quarks) , as illustrated in Fig. In this case, only the lower branch is occupied, and the total energy 
lowers for formation of the density wave with wavenumber 2k F due to the nesting effect [l2l.ll3| . In the uniform 
3-D system we discuss in this article, however, the wavenumber dynamically depends on the balance between the 
kinetic and interaction energies, and becomes smaller than 2k F '- the spatial modulation due to DCDW makes 
the kinetic-energy loss, and the energy gain is generated by the deformation of the single-particle spectrum 
which originates from the p-h interaction. This situation is essentially the same as SDW in 3-D electron system 
discussed by Overhauser [lj|. In the next section we demonstrate the actual manifestation of DCDW by taking 
a definite model. 



±^+M 2 +?/2+Ei(p) \ 

Pic+iPy , (8) 

— ^==^3 A ±} (p + q/2). 0) 



III. APPLICATION TO THE NJL MODEL 



Since DCDW defined by eq. is associated with the chiral condensation, we consider the moderate density 
region where nonperturbative phenomena are expected to remain even in quark matter. We here employ the 
NJL Lagrangian with Nf — 2 flavors and N c = 3 colors 0, 0] to describe such a situation, 

Cnjl = $(ip - m c )V> + G[ W) 2 + (^75T^) 2 ], (11) 

where r is isospin matrix, and m c is the current mass, m c ~ 5MeV. 
We assume the mean-fields in the direct (Hartree) channels, 

($1/)) = Acos(q-r) 
(^757-3^) = Asin(q-r), (12) 

where we fix the isospin direction to T3 because of degeneracy on the isospin hyper sphere; the other mean-fields 
vanish consistently, (ipi^fsTiip) = (^175X2*0) = 0. In other words, it is assumed that DCDW is a charge eigen 
state, and there is no amplitude which mixes states with different charges [2£|. As for the Fock exchange terms, 
we have briefly examined them in Appendix [51 and shown that the tensor exchange terms might affect DCDW 
implicitly through the deformation of one-particle spectrum. In the present study, however, we will treat only 
the direct terms since the exchange terms correspond to pure quantum processes and have less contribution than 
the direct ones. It is interesting that the configuration i|12[) is similar to the pion condensation in high-density 
nuclear matter within the a model, suggested by Dautry and Nyman poj. where a and tt° meson condensates 
take the same form as eq. <|12[) . It may implies a kind of quark- hadron continuity |2l|. 
Within the mean-field approximation, the effective Lagrangian becomes 

M 2 

C MF =tjj[ip+ pi 7o - M{cos(q • r) + 17573 sin(q • r)}] ip - — , (13) 
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(a) Pz (b) Pz (c) J> 




FIG. 2: Fermi surfaces of the spectrum for q//z. pt = \/pl + Py- The outer (inner) closing line corresponds to 
E-(p) (E+(p)). (a) for M > q/2; the minimum of E-(p) is at the origin, (b) for M < q/2; there are two minima at 

points (pt,Pz) = (0, ±-\/ (?/2) 2 — M 2 ) . (c) for M < q/2 and /i = q/2 + M; the minor Fermi sea vanishes. 



where we have introduced the chemical potential /i, and taken the chiral limit (m c — 0) assuming /i >> to c . 
Since only the difference between u- and d- quark is the sign of the wavenumber vector (q <-> — q) due to the 
isospin matrix T3, the single-quark spectrum takes the same form as in eq. (J5J). Thus we need not distinguish 
two flavors in the energy spectrum, though the eigen spinors depend on the sign of q. 

Hereafter, we take the direction of the wavenumber vector parallel to the z axis, q = (0, 0, q) without loss of 
generality, and show the Fermi surface for various values of /x, M, and q in Fig. [21 

In the density-wave state, all the energy levels below the chemical potential are occupied in the deformed 
spectrum. Accordingly the thermodynamic potential density at zero temperature becomes 

f d 3 77 M 2 

"tot = N f N c / 7^4 E li E '(P) " M> * 0* - E s (p)) E S { P )] + 4g 

S — ± 

M 2 

= flfer + n v! » c + — , (14) 

where f2f or (^vac) denotes the Fermi-sea (Dirac-sea) contribution. We can see that the finite wavenumber effect 
enters only through the deformation of the energy spectrum, and gives non-trivial contributions to the behaviour 
of the chiral condensation or the dynamical mass M. The energy gap between the two branches is generated 
by the dynamical mass which comes mainly from the Dirac sea, and the energy gain due to the density wave 
(to a finite q) comes essentially from the Fermi sea which is responsible for the finite baryon-number density, 
thus the DCDW is produced cooperatively by the Dirac and Fermi seas. 

Since the NJL model is unrenormalizable, we need some regularization procedure to evaluate the negative- 
energy contribution Q vac . Because of the spectrum anisotropy we cannot apply the momentum cut-off regular- 
ization scheme. Instead we adopt the proper-time regularization (PTR) scheme We show the result (the 
derivation is detailed in Appendix lA"|l . 



_ 7 / dr / dp 

* 'va.i 



8^/2 y i/A2 r 5/2 J ^ 2n 



" Slrei, (15) 



where A is the cut-off parameter, and we subtracted an irrelevant constant O re f in the derivation. All the 
physical quantities should be taken to be smaller than the scale A in the following calculations. 



A. Phase transition at zero temperature 



To investigate threshold density for formation of the density wave at T = 0, we expand the potential (|14|l up 
to the second order in q, and examine the sign of its coefficient, 

fi tot = ttL + ^(/3fcr+/3vac)<Z 2 + 0(<Z 4 ) (16) 
B 2 VLr M 2 

/3 fer ee ^H^o = -NfN^H^/M) (17) 
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/3vac = -j^U^o =N f N c — J(A/ 2 /A 2 ), (18) 



where J(x) = x drexp(— t)/t, and -ff(x) = ln(a; + \/# 2 — !)• The coefBcient of the second-order term in 
fifer is always positive for finite dynamical mass M ^ 0, while the counterpart of f2 V ac is negative, indicating 
that the Dirac sea is stiff against the formation of the density wave. In contrast, the Fermi sea favors it as 
mentioned in the previous section. The total coefficient, /3 tot = /3f er + Pv&c, depends on the dynamical mass and 
the chemical potential for fixed A, as shown in Fig.|3 For larger values of the chemical potential in Fig.^b), 

(a) (b) 



(i/A=0.5 




0.00 0.05 0.10 0.15 0.0 0.2 0.4 0.6 

M/A M/A 



FIG. 3: The dynamical-mass M dependence of the coefficient f3 to t for various values of the chemical potential, (a) for 
At/ A = 0.37 - 0.39, (b) for p/A = 0.5 - 0.7. 

/?tot becomes negative and reaches its maximum at a finite M. As for the small chemical potential p/A < 0.38 in 
Fig. EH a), 0tot never become negative for any value of M. It leads to a rough estimation of the critical coupling 
constant, GA 2 ~ 4.63, which is the value to occur the usual chiral condensation (q = 0) at fi/A = 0.38 in the 
PTR scheme. 

The magnitudes of M and q are obtained from the minimum of the potential (|I4|I at T = 0, and their values 
satisfy the stationary conditions, dtttot/dM = dfltot/dq = 0. 

Fig. 0] shows contours of S! to t in M-q plane for a given chemical potential, where the parameters are set as 
GA 2 = 6 and A = 850 MeV, which are not far from those for the vacuum (p = 0) [l^. The crossed points 
denote the absolute minima. There appear two critical chemical potentials p — p c i,p C 2- for the lower densities 
(Fig. |2£a)- (b)) the absolute minimum lies at the point (M ^ 0, q = 0) indicating a finite chiral condensation. At 
P = A*ci (Fig-EI c )) the potential has the two absolute minima at (M 0,q — 0) and (M y^0,q^ 0), showing 
the first-order transition to the DCDW phase, which is stable for fi c \ < fi < p C 2 (Fig. @Jd)-(e)). At fi = p c % 
(Fig. 01f)) any point on the line M = and a point (M ^ 0, q ^ 0) become minimum, and thereby the system 
undergoes the first-order transition to the chiral-symmetric phase which is stable for fi > /i C 2. 

The Fig. [3] shows the behaviors of order-parameters M and q as functions of /j, at T = 0, where that of M 
without the density wave is also shown for comparison. It is found from the figure that the magnitude of q 
becomes finite just before the critical point of the usual chiral transition, and DCDW survives at the finite range 
of n (/i c i < fx < /U c 2) where the dynamical mass is reduced in comparison with that before the transition, and 
decreases with fi. On the other hand, the wavenumber q increases with /x, but its value is smaller than twice the 
Fermi momentum 2/cp(~ 2[i for free quarks) due to the higher dimensional effect; the nesting of Fermi surfaces 
is incomplete in the present 3-D system. Actually, the ratio of the wavenumber and the Fermi momentum (at 
normal phase q = M = 0) becomes q/kp = 1.17 — 1.47 for the baryon-number densities pb/po = 3.62 — 5.30 
where DCDW develops. The baryon-number density is shown in Fig. as a function of p for the normal and the 
density wave cases. The jumps of the baryon-number density reflects the first-order transition. In the DCDW 
phase, the relation q/2 > M is retained, and the Fermi surface looks like Fig. Etb) . 

Here we show the coupling-strength dependence of the critical chemical potentials p c i,c2 in Fig. d including 
the semi-empirical value, GA 2 = 6.35 (A = 660.37 MeV), to reproduce the pion decay constant / ff = 93 MeV 
and the constituent-quark mass ~ 330 MeV, one third of the nucleon mass in the vacuum. The range of DCDW 
phase between p c \ and p C 2 starts to open at GA 2 ~ 4.72, and broadens with increase of the coupling strength. 
It should be noted that the effective potential HI4[) can be scaled by the cut-off A and thus the dimensionless 
coupling GA 2 becomes only one parameter to determine whether or not the phase transition itself occurs in the 
present model. 
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M/A M/A M/A 

FIG. 4: Contours of fltot at T = are shown in M-q plane as the chemical potential increases, (a) — > (f). The cross in 
each figure denotes the absolute minimum. 



Order Parameters at 1=0 
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FIG. 5: Wavenumber q and dynamical mass M are plotted as functions of the chemical potential at T = 0. Solid (dotted) 
line for M with (without) the density wave, and dashed line for q. 



It should be kept in mind that the order of the transitions may depend on the parameter choice of G and A, 
and may also on the regularization scheme. 



B. Magnetic properties 

Using the eigen spinors in eqs. (fH|>- (|10|) . we can calculate various expectation values with respect to the DCDW 
state. For a operator, O, which does not depend on the spatial coordinate, its expectation value becomes a 
simple form: 

JUL $t (p) e -375q-r/2 0e - J r 375 q-r/2 ^/ ^ 
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FIG. 6: Baryon number density as a function of fi. po = 0.16fm 3 : the normal nuclear density. Solid (dashed) line is for 
the finite q (q — 0). 
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FIG. 7: Critical chemical potentials fi c i,c2 are plotted as functions of the dimensionless coupling GA 2 at T = 0. DCDW 
phase appears between fi c i and /i C 2- The critical coupling at /x c i = /j, C 2 is estimated to be GA 2 ~ 4.72. 



We can confirm that baryon-number density, O = 1, is still constant even in the density wave state: summation 
of quasi-particle state in momentum space, 



PB 



d 3 p 

(27T) 5 



(tp ip)ip{p)) = constant. 



On the other hand, the spin expectation value, O — 707573/2 = S z , vanishes in each flavor, 

1 ftf^M, 



because the stationary condition for the wavenumber q is proportional to the expectation value: 

= ¥q «<E,). 



(20) 



(21) 



(22) 



Here we show an interesting feature of DCDW: a spatial modulation of the anomalous magnetic moment. 
The Gordon decomposition of the gauge coupling term gives the magnetic interaction with external field f" 
in the form, <?i(e* /2M)(i/j<7 ^v^F^ ', where is a form factor and e* an effective electric charge. The operator 
of the magnetic moment for the z component is defined by O = 70^12, which is not commuted to 75, 



e ir 37 5qT/2 7o(ri2e -ir 376 q.r/2 = CQs(q . ^ _ gin(q . ^ 



(23) 
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and then its expectation value is given by 



(ip(r)<rui/)(r)) = (70CT12) cos(q • r) - £{73} sin(q • r), (24) 
d 3 p 2M 

(73) = 0. (26) 



where (70^12) = / 7^— 3 — 7==== [n + (p) - n_(p)] , (25) 



The function n± (p) is the momentum distribution for the eigen state corresponding to E± (p) . The expectation 
value (70C12) is proportional to an asymmetry of the momentum distribution in n±(p). In DCDW phase, 
the asymmetry becomes finite as shown in Fig. |2 and thus the magnetic moment is spatially modulated 
with wavenumber q. The expectation value (73) vanishes in each eigen spinor: ul^aii-t = v±73i>± = 0. We 
also confirmed that the other components of the magnetic moment, 700-23,13, vanished analytically after the 
integration in momentum space. 

The equation H24|) shows that the amplitude of the modulated magnetic moment depends on the dynamical 
mass, reflecting delay of chiral restoration due to the presence of DCDW. The magnetic order of DCDW should 
have some observable consequence of compact stars with quark cores. We estimate the magnitude of the 
amplitude (|2*5|l at ps/po = 3 ~ 4: the tensor expectation value per quark is calculated to be (7oci2)/(7o) = 
0.1 ~ 0.3. Thus a local magnetic flux $ induced by DCDW, 

-(H) ^7>- <*> 

amounts to O(10 16 ) Gauss, which is comparable with observed values in magnetars [2(|. The flux strength on 
the star surface from the quark core might be smaller, since it is given by summation of the quark-magnetic 
moment (12411 modulated rapidly with the wave length g/27r ~ O(10) fm, nevertheless, contributions from near 
the quark-core surface may remain without the cancellation. 



C. Correlation functions 



In this section, we consider scalar- and pseudoscalar-correlation functions, II SiSp (fc), at the chirally restored 
phase, and discuss their relation with DCDW at T = 0. The correlation functions depend on an external four- 
momentum k = (fco,k), chemical potential, and also the effective quark mass for a given chemical potential. 
In the static limit fco — * 0, the correlation functions have a physical correspondence to the static susceptibility 
functions for the spin- or charge-density wave [TT| . while they have no primal singularity at |k| = 2k p reflecting 
the higher dimensionality. Note that the functions have a differential singularity at k — 2kp. 

In general, a pole of the correlation function implies a second-order phase transition. As for the case of 
DCDW, the effective potential analysis in the previous section shows the first-order transtion. Nevertheless, as 
shown in Fig.QJf), the potential barrier from DCDW to the normal phase at the critical point pL = \i c % is very 
small, and it is expected that the correlation function has some imformations of the transition point. 

We evaluate effective interactions, r s>sp (fc), within the random phase approximation |17lll9( . which are related 
to the correlation functions, i.e., 2GII SjSp (fc) = r S)Sp (fc)II^ sp (fc): 



*r S)PS (fc) — 



2Gi 



2Gn s % s (/c)' 



(28) 



where 11° DS (fc) is the polarization function in medium, which has form at the static and chiral limits (see 
Appendix tO): 



n°(|k|) = n°(|k|) = 



N f N c 



At: 2 
N f N c \k 
4tt 2 



(A 2 - 2k F ) - 2Ar / A c ik 2 /(k 2 )| 



M-*0 



k F log 



(2ki 



\2k F - |k| 



M log |(^. 



|k|)(2fcf-|k|) 
Ikl 2 



(29) 



Inverse of the correlation function in the chiral limit corresponds to the coefficient of M 2 of the effective potential, 

1. 



fitot = O tot | M ^ + -r- s L (q)\ M ^ a M 2 + 0(Af 4 ). 



(30) 
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FIG. 8: Function, l/r ps (|k|), is plotted for various values of fi, fi/A = 0.4911, 0.5, 0.53254(= (la), and 0.55, for m c = 
(a) and m c — 10 MeV (b). In the case of fi/A > 0.4911, the mass-gap equation has a extremum solution M = 0. 



From behaviors of the function r ps (|k|) _1 shown in Fig.[S{a), it is found that the function takes the lowest value 
at a finite external momentum (|k| ~ l.bkp), and thus a finite wavenumber q gives the lower potential energy 
in eq. (|3U|) . Assuming the second-order transition, the critical density fit and wavenumber |k t | are determined 
from a following simultaneous equation: 

r sp (|k|; M )- 1 =0 ) and ^^^1 =0. (31) 

A numerical calculation in the chiral limit gives fit{= kp) = 0.5320A and |k t | = 1.498fcp, which almost coincide 
with the correct result fi c2 = 0.53254A and q = 1.469fc F (where k F = y/ fi 2 c2 ~ M2 '-> M = 0.034A) from the 
effective potential analysis. A complex structure of the effective potential makes a little difference between |k t | 
and q at fi = /i C 2- 

The above argument might also be available even for the case of a finite current-quark mass, m c ~ lOMeV: 
the effective potential for a small current-quark mass (m c << fi) is approximated to 

n tot ~ n tot \ M -> mc + h:^{ q )\ M ^ mc {M - m c f + o[(M - m c ) 4 ]. (32) 

The Fig. EJb) shows that behaviors of the coefficient function has little shift from that of the chiral limit, and 
thus suggests a DCDW transition with a small current-quark mass as in the chiral limit. 



D. Phase diagram on the \i-T plane 



To complete a phase diagram we derive the thermodynamic potential at finite temperature in the Matsubara 
formalism. The partition function for the mean-field Hamiltonian is given by 



Z. 







DipDip exp / dr / d r < ip 



id + M exp (Hsq • r) - 70/i 



(33) 



where /? = 1/T, and d = — 7o<9 r -H7V. Taking the Fourier transform of the spinor with the Matsubara frequency 
u> n , the partition function becomes 



^(r,-fr) = e -^75q-r/2 T ^ J 



d 3 k 

(2^)3 



■0(k, n) 



J] {(^ n + M ) 2 -£; s 2 (k)}^xexp 



k,n,s— ± 



f-2 



M 
4G 



V/3 



(34) 
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FIG. 9: The wavenumber q and the dynamical mass M are plotted as a function of /i at finite temperature. The solid 
(dashed) line corresponds to q (M) at T/A = 0.02 (thick lines), 0.03 (thin lines). 




FIG. 10: A phase diagram obtained from the thermodynamic potential eq. I|M5^ . The solid (dashed) line shows the 
chiral-transition line without (with) the DCDW. The shaded area shows the DCDW phase. 



where V is volume of the system. Thus the thermodynamic potential flp per unit volume is obtained, 



Sl fi (q,M) = -T log Zp(q,M)/V 



e -/3(E s (k)-p) + x 



+ 1 



E, 



1 M 2 
(k) } + 4G (35) 



where we have utilized a contour-integral technique for the frequency sum to get the final form. 

From the absolute minimum of the thermodynamic potential (|35|l . it is found that the order parameters at 
T ^ behave similarly to those at T = as a function of fi, while the chemical-potential range of DCDW, 
fJ-ci{T) < fi < ^ c2 (T), gets smaller as T increases. The Fig.|5]sliows the order parameters at finite temperature. 
The discontinuities of the order parameters reflect the two absolute minima at the critical chemical potential 
Hci,2, and it indicates a first-order transition. Thus the region of DCDW in the \x-T phase diagram is surrounded 
by the first-order transition lines. 

We show the resultant phase diagram in Fig. ^3 where the usual chiral-transition line is also given for 
reference. Comparing phase diagrams with and without q, we find that the DCDW phase emerges in the area 
(shaded area in Fig. lll)|) which lies just outside the boundary of the ordinary chiral transition. We thus conclude 
that DCDW is induced by finite-density contributions, and has the effect to extend the chiral-condensed phase 
(M ^ 0) to a low temperature (T c ~ 50MeV) and high density region. The above results suggests that QCD at 
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finite density involves rich and nontrivial phase structures, as well as color superconducting phases. 



IV. SUMMARY AND OUTLOOKS 



We have discussed a possibility of the dual chiral-density wave in moderate-density quark matter within the 
mean-field approximation, employing 2-flavor and 3-color NJL model. The mechanism of the density wave is 
quite similar to the spin density wave in 3-D electron systems; the total-energy gain comes from the Fermi-sea 
contribution in the deformed spectrum, while its amplitude has the different origin corresponding to the chiral 
condensation from the Dirac-sea effect. 

In this paper, we have considered only the direct channels (Hartree terms) of the interaction. If the exchange 
channels (Fock term) are involved, there appear additional interaction channels by way of the Fierz transfor- 
mation 01 • I n particular, self-energies in axial- vector and tensor channels related to a ferromagnetism |rj l25j 
might affect the density wave through nontrivial correlations among them. Interactions in the p-p channels are 
also obtained by the Fierz transformation, and their strength is smaller than that of the direct channels by the 
factor of 0(1/N C ). Because the Cooper instability is independent on the strength of the interaction, it is inter- 
esting to investigate the interplay among the density wave, superconductivity p. IsL |27| . and the other ordered 
phases, e.g., chiral-density waves mixing isospins which may cause a charge-density wave due to difference of 
electric charges of u- and d-quarks, as future studies. 

Finally, it is worth mentioning about correlation functions (fluctuation modes) on the density-wave phase, 
which give the excitation spectrum, and are important for the dynamical description of the phase. In par- 
ticular, Nambu-Goldstone modes are essential degrees of freedom for low-energy phenomena, and may bring 
some observable consequences, e.g., slowing down of star cooling through enhancement of specific heat due to 
fluctuations of such low-energy modes. 
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APPENDIX A: REGULARIZATION OF Q 



We regularize the Dirac-sea contributions to the potential, f2 V ac, by applying the schwinger's proper-time 
method. f2 V ac can be described in the form of the one-loop order contribution, 

Ail = n vav (q, M) - fi N = - / -£jL £ log £l (Al) 

Jc H 2 *") 4 fr± d n 

with D ± =kl-El(k) and D N = fc 2 , - k 2 - to 2 , (A2) 
where Cl^ is the normal vacuum contribution. Using the identity for G G R 

/>oo 

(G + iri)- 1 = -i ds e i( - G+ir >> , (A3) 
Jo 

we find 

D± + in f°° ds ( a i( D±+iv ) s _ j(Dn+iv)s 



log^^- - = -/ — ( e K"±+iV)s . (A4) 

D N + ir) J s V / 
By way of the Wick rotation which is done simultaneously for ko integration of fi rac and Qnj 



c *(2^) 4 Jo s 

/oo i4 f, f 'Oo i 
^ / H£ expi{-fc| - kl - {^WTJP±a/2f} S 
-oo (2tt) 4 Jo s I 
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expi{— k E — k\ — k z — m 2 }s] 



50 dr 

(27T) 4 J T 



dk z 



exp{-fc| - k 2 t - (y/kJ+m±q/2) 2 }T 
- exp{-fc| - k 2 t - k 2 z - m 2 }r] 



dr 



8tt 3 / 2 J 2tt 7 r 5 /2 



exp{-( v //j2 + M 2 ± g/2) 2 r} - exp{-(fc 2 + m 2 )r} 



(A5) 



The above integration of r has singular at r = and thus not well defined. The proper-time regularization 
is to replace the lower limit of r by the cutoff 1/A 2 , 



dr 



dr, 



(A6) 



1/A 2 



the A corresponds to a momentum cutoff. 

Eventually we obtain the regularized potential from the Dirac sea, 



AO 



8tt 3 / 2 



dfc, 
2vr 



dr 



1/A 2 T 



5/2 



cxp{-(v/fc 2 + M 2 + <?/2) 2 r} + exp{-( v /fcfTM 2 - <z/2) 2 r} 
-2exp{-(fc 2 + m 2 )r}] . 



(A7) 



The normal vacuum contribution Q n has an explicit form, 



n N = N r N 



1 



dr 



8tt 2 ./ 1/A 2 r 3 



A' 



16tt 



1/A 

*{<- 



■m 2 + l)e-™ 2 



4 r(o,77V)} 



(A8) 



where in = m/A and F(a, z) = J^ 00 drr° 1 exp(— r) is the incomplete gamma function. 



APPENDIX B: FOCK EXCHANGE CONTRIBUTIONS 



We briefly examine how the Fock exchange terms of the NJL interaction affect DCDW. After the Fierz 
transformation, one can find the exchange interaction terms, discarding color-octet contributions 19] : 



T 



1 

8Ac 



2 (^) 2 + 2 (tpi-f 5 T^) 2 - 2 (^tiP) 2 - 2 (^75 V) 2 



(Bl) 



The overall factor 1/8N C indicates that the exchange terms are less relevant in comparison with the direct terms. 

The first two terms come to be added to the Hartree terms, and contribute to DCDW through changing the 
effective interaction in these channels by a factor. 

The third and fourth terms have opposite sign of interactions, and thereby they can not gain condensation 
energies. 

The fifth one has an expectation value in its temporal term, which corresponds to the baryon number density, 
(i/j^i/j) ^ 0, and is renormalized into the chemical potential |23,[24|. Thus it might not change DCDW quali- 
tatively, except for changing value of a bare chemical potential for a fixed baryon number density. The other 
spatial terms vanish self-consistently in the present formalism 0j . 

The sixth term, the axial-vector interaction, seems to contribute DCDW because the Weinberg transformation 
leaves it unchanged, but scalar and pseudo-scalar terms for DCDW changed to an another axial- vector mean-field 
with isospin T3. However, it is found that this term does not affect DCDW by itself: the Hartree- Fock free- 
energy with both the axial-vector mean-field Va and DCDW (the wavenumber q) is given, after the Weinberg 
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transformation, by 
Hhf = 



^ [p 2 



M 2 



7375 (Va + 739/2)] ip- 



M 2 



4G S 4G A 



ip u [p 2 + M 2 + 7375^] tp u + / i>d [p 2 + M 2 + 7375S] Vrf + 77^ + 



M 2 (A + B) 



4G S 16G A 



(B2) 



A = + g/2, and B = V A - g/2, 



where ^ u ,d are spinors of u,d-quark, and Gs,a are reduced coupling constants in scalar and axial- vector channels. 
The equation l|B2(l shows that the minimum of the effective potential is always given by A + B — 2Va — 
because the first two terms have the same structure in their energy spectrum with respect to A, B independently 
of the sign of A, B, and implies, therefore, that the expectation value of the axial-vector channel vanishes. 

The last two terms, the tensor interactions, might give a non trivial contribution to DCDW through an 
feedback effect of the nonzero expectation value of the tensor operator, (700-12), see eq. (J2U, although their 
interaction strengths are small by the factor 1/8N C . This problem is left for a future study. 



APPENDIX C: SCALAR AND PSEUDO-SCALAR SCATTERING AMPLITUDES 



Following Nambu |17| , we consider the quark-quark scattering matrix generated by the chain diagram in the 
pseudo-scalar channel. Then the polarization function, — iIIp S , is given by 



with the quark propagator in medium, 

Sip) = 



1 7T 

+ i — {i> + m)d{n - E p )6{ Po - E p ) 



f — m 



{j> + m)S(p) = {j> + m) S F (p) + S D (p) 



Then the scattering matrix M?, can be written in the form, 



zM 3 p 3 s (fc 2 ) = ( i75 )T 3 



There are three kinds of contributions to Ilp S (fc 2 



2iG 



1 - 2GH° (fc 2 ) 



(«75)7 3 . 



e° (fc 2 ) = n^(fc 2 ) + n^(fc 2 ) + 2irr (fc 2 ) 



with 



in^ s {k 2 ) = -w f N c 



First, we consider the vacuum contribution, 



zn£ F (fc 2 ) = 2N f N c 



d 4 p 



(2n) 



(p + \k) 2 — m 2 (p — ifc) 2 — m 2 



-2N f N c k 2 



d 4 p 



(2tt) 4 [(p+ ^k) 2 -m 2 ][{p~ \k) 2 ~ 



The integrals in the first term is easily evaluated to get 

d 4 p 1 

(2tt) 4 ( p ± \kf-m? ~ 



/ 



(CI) 



(C2) 



(C3) 



(G4) 



(C5) 



(C6) 



(C7) 
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in the proper-time representation. The integral in the second term is denoted by I(k 2 ), 

d 4 p 1 



Using the Feynman's trick, 



Jo J (2tt 
and introducing the proper-time r, we find 



(2tt) 4 [(p + ifc)2 - m 2 ] [(p - ifc) 2 - m 2 ] ' 

d 4 p 1 

(2tt) 4 [p 2 - to 2 + (k 2 + 2pk)t} 2 
d A p d 1 



dt 



The vacuum contribution is summarized as follows: 



iNfNc 2 
~^ 2 ~ k 



f 1 , /*°° dr 
/ dt / — ( 



-(m 2 -fc 2 t(l-t))7 



Secondly, let us consider U^ s F (k 2 ), 



(ko + E p _^) -(p + ^-m 2 ^ 

Taking the static limit k — > 0, we have 

. f d 3 p pk 1 

"J w, 



-0(M-£p-ik)- 



-ffl"/(|k|) - 2JV>JV C * 



2A^/7V c i / da; 



r) 3 2p • k + k 2 Ej 
p 2 dp 



-6(j* - E p ) 



2px+ |k| 



E, 



■6(n -E p ) 



.NfN c 



+ 



Air 2 

N f N c \k 
8?r 



PF[i. — m In 



/i - m + 



m 



| k |2xV 2 ' 



pdpln 



2p+|k| 



2p-|k| 



E„ 



(C8) 



(C9) 



(CIO) 



(Cll) 



(C12) 



(C13) 



The integral over the Fermi sea can be analytically performed, but we do not write it down here because of its 
complexity. 



Finally, we calculate H ps D (k 2 ), 



-m° D (k 2 ) = - 



d 4 P 

(2^) 4 
d 3 p 

¥ 



J (27T 



4iV>iV c [-p' + m z + -k z )S D [p+-k\S D [p--k 



ANfN c 



-ko + E p ) + p + -k + m z + -k 



9(ji - E p+li )5(k + E P - E p+k )9(fi - E p ). 



In the static limit k a — > 0, 



^ (|k| 



M- | m + — 



k 2xV 2 



(C14) 



(C15) 
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which exactly cancels the imaginary part arising from IlSf (|k|). 

Collecting them together, we have the denominator of the scattering amplitude in the static limit, 

i-2Gn° s (|k|) = i-2G(n^(|k|) + 2n^(|k|) + n^(|k|)) 

,N f N c fdr, T 



1 - 2G- 



4tt 2 



AGiN f N c k 2 I(k 2 ) 



G 



N f N c 



Pf^ — rn 2 In 



(pf + v\ |k| f p 

\ m 



pdplrt 



2p 



2p-\k\ 



En 



(C16) 



On the other hand, the gap equation in this case reads 



m = m c + 



1 

2^ 



:GN f N c m 



dr 



G 



N f N c 



PF[i — rn 



2 In 



Pf + M 



(C17) 



Thus, we find at a stationary point (at a solution of the gap equation) 

i-2Gn° s (|k|) 



rn 



4GiN f N c k 2 I(k 2 ) - G^-^\k\ I pdp\n 



2vr 2 



2p 



2p - |k| 



E,, 



(C18) 



In the similar way, the scattering amplitude of scalar channel is given by 

l-2GII°(|k|) = 1-2G^^ / %e~ m2r ~4GiN f N c (k 2 +4m 2 )I{k 2 ) 

47T J T 



G 



N f N c 



2 (PF+V 

Pf^ — rn in 



V m 



k 2 + 2m 2 
2|k| 



pdpln 



2p 



2p- 



2GN f N c i 



m 
2vr|k| 



H - ^m 2 +k 2 /A\ 9 {2k F - |k 



(C19) 
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